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Abstract 

Mean field and exact diagonalization studies on the quasiparticle excitation 
of an SO(5) symmetric two-leg ladder system are reported. It is shown that 
the energy gap in the quasiparticle excitation is caused by the formation of 
rung singlet states. We find that shadow peaks can occur above the Fermi 
surface with antiferromagnetic electron correlations involving only rungs in 
the spin ladder. 
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The cuprate materials of high critical superconducting temperatures exhibit antiferro- 
magnetism near half-filling and superconductivity away from half filling. Zhang [|T],3 sug- 
gested that the two phenomena are manifestations of the same thing and can be explained 
in a unified framework based on SO (5) symmetry. Recently Scalapino et al. || presented a 
simple 5*0(5) model on a two- leg ladder system which can be served as a toy model to inves- 
tigate various physics involving the SO (5) symmetry. In this paper we discuss quasiparticle 
excitations for the SO (5) two-leg spin ladder from both analytic and numerical results. We 
derive a mean field quasiparticle energy dispersion relation. For a through comparison one- 
particle spectral function is calculated by applying a Lanczos exact diagonalization method 
to a ladder of 2 x 6 sites. We find that the energy gap in the quasiparticle excitations is 
largely contributed by the formation of rung singlet states. We also find the appearance of 
shadow peaks with Heisenberg interaction along rungs alone. 

The 5*0 (5) symmetric two- leg ladder Hamiltonian is given by [[3] 

H = ~t\\ Y,( c l\* c i+i\* + H - c -) - t± ^( c Lc i2CT + H.c.) 

i\cr ia 

+UJ2 (that " I) (n iH - ~) + VX(rm - l)(n i2 - 1) 

+ ■ S i2 ~ y-^riiXa (1) 

i i\cr 

Here c\ Xa creates an electron with spin a on the z-th rung of the A-th leg with i = 1, . . . , L 
and A = 1,2. n is the number operator, n iAcr = c\ Xa c iXcr , rii\ = Y^a n i\<n an d S is the spin 
operator, = \c\\ a a 'ap c i\p- t\\ is the hopping integral in the leg direction; t±, the hopping 
integral in the rung direction; U, the on-site Coulomb interaction; V, the near-neighbor 
Coulomb interaction on a rung; J, the Heisenberg exchange interaction on a rung; and fi, 
the chemical potential. The Heisenberg interaction of the SO (5) ladder Hamiltonian allows 
both singly and doubly occupied sites, while that of the t-J ladder Hamiltonian allows only 
singly occupied sites. The constraint for the interaction strengths, J = 4(U + V) in Eq. ([!]) 
above is required for the 50(5) symmetry. U and V can be repulsive or attractive potentials. 
In the present study we will consider only the case of repulsive potentials, U > and V > 0, 
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which is in the region of the Eq spin gap phase (represented by the tensor product of rung 
singlet states in the strong coupling limit) defined by Scalapino et al. || 

The weak coupling (that is, U, V <C tii ~ t±) may allow linearization of the Coulomb 
repulsion terms 

1\ / 1 
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and 

(t^i - l)(7i i2 - 1) 

= (riii - l)(n.j 2 - 1) + (n i2 - l)(n« - 1) - (n%i - l)(n i2 - 1) 

Taking into account (n^t) = (^aj.) = ^ with 5, the doping rate, which neglects hole 
density fluctuations, we write the Coulomb repulsions (n^T ~~ |)( n iA4. — §) = ~ | n »A + f — T 
and (riji — l)(7ij 2 — 1) = — 5(nji + rij 2 ) + 25 — <5 2 , and add them to the chemical potential 
term in Eq. ([[]). The Heisenberg interaction can be written [|J as 

3 3 
S,i • S i2 = --[x*i 2 ( c liT c i2 T + c !i| c i2|) + H - c -] + g^il 
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g [^il2 ( C ilt C i2i — C ill C i2l) + H.C.] 

where the hopping order parameter is Xn2 = ( c lit c i2t + c lu c i2j) an d the singlet pair order 
parameter, A il2 = (% T c i2i - c iU c i2T ). We take Xa 2 = X and A a2 = A by neglecting spatial 
fluctuations of both the amplitude and the phase. The mean field Hamiltonian is then in 
momentum space, 

H = E(~ 2t ll cos k * ~ v)A\c c k\o ~ t± E(4i CT c fc2(J + H.c.) 

k\cr kcr 
3 J i , 

jT~ 5Z[^*( c fciT c fe2T + c fcii c fc2j.) + ^*( c fciT c -fc2| — c kii c -k2^) + H.c] (2) 
k 

The quasiparticle energy dispersion is readily obtained from Eq. (0) above, 
E k = ± 



C — 2*11 cos ^ - //) ± (t_L + -^J 
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The mean field order parameters Xi A> and the chemical potential \x are obtained from 
self-consistent equations. For the time being we take the values of %, A, and fi obtained 
from our Lanczos calculations on the 2x6 ladder with two doped holes. The quasiparticle 
energy dispersion is shown in Fig. [I]; the bonding band (denoted as B) has its minimum 
energy at k x = with the Fermi surface at kp ~ while the antibonding band (denoted as 
A) has its minimum energy at k x = with the Fermi surface at kp ~ |. Splitting between 
the two bands, that is, the removal of degeneracy is caused by hopping (overlap) integral 
t± in the rung direction. The bonding band is pushed down with more occupied electrons, 
while the antibonding band is pushed up with less occupied electrons. We note that with 
kp > kp the Luttinger sum rule is satisfied, that is, kp + kp = (1 — 5)n. The dashed parts 
of the dispersion curves in Fig. [I] represent the shadow bands in our two-leg ladder system 
which is similar to the ones discussed in two-dimensional planar systems As shown 

in Eq. (^) the shadow band in Fig. [I] is largely contributed by the Heisenberg interaction 
in the rung direction, while the shadow band in the two-dimensional planar systems 
is caused by antiferromagnetic correlations whose strengths are equal in both the horizontal 
and vertical directions. 

To examine the single particle excitations, we write in the case of two-hole doped system 
of 2 x L sites @, 

A e {^u) = £ l^lcLl^ 2 )! 2 5 ^ - E " L " + E o L ' 2 + A*) ^ 

a 

for the particle spectral function {uo > 0) and 

A h (k,u) = £ \(K L - 3 \c^l L - 2 )\ 2 - £f' 3 + Ef~- - „) (4.6) 

a 

for the hole spectral function {u < 0). Here ) is the a-th eigenstate in the subspace of 
iV electrons with the eigenenergy ■ The chemical potential is defined as [i = ^(E$ L ~ l — 
Eq L ~ 3 ). We calculate the spectral function by using both the Lanczos method and the 
continued fraction approach 0. The results are shown in Fig. ^|for L — 6, t\\ — t± — 1 and 
V = 0, as a function of U and J. The coupling strengths, U and J are related to each other 
and J increases faster than U due to the SO (5) constraint, J = 4(C7 + V) = AU . 



We now investigate the spectral functions by varying the values of U (and consequently 
J). First the computed free particle spectral functions for the case of U = J = are shown 
in Fig. 0(a). The Fermi momentum of the bonding (k y = 0) band is kp = and that of 
the antibonding (k y = tt) band is kp — |. Below the Fermi surface, we observe sharp and 
well-defined quasiparticle peaks for both the bonding and the antibonding bands. All the 
single particle states below the Fermi surface are occupied, and there exist no spectral peaks 
above the Fermi surface. The computed spectral functions shown in Fig. ^|(a) are seen to 
agree well with the free particle energy dispersion, e k = —2t\\ cosk x ± t± in Eq. (|]). 

The spectral functions for the case of small coupling strengths, U — 0.1 and J = 0.4 are 
displayed in Fig. |2|(b). The predicted Fermi momenta are found at kp ~ ^ and kp ~ | 
in agreement with the mean field result given by Eq. (^|). The calculated positions of the 
quasiparticle peaks below the Fermi surface are also in good agreement with the analytic 
mean field result of Eq. (|3]). Although the interaction strengths are small, interestingly 
enough there appear spectral peaks above the Fermi momentum (surface). Such peaks are 



the shadow peaks which were also observed by others [ID|. The peak in the bonding band 
with momentum (tt, 0) above the Fermi surface is a shadow of the peak belonging to the 
antibonding band with momentum (0, tt) which is below the Fermi surface. The peak with 
momentum (y-,7r) is a shadow of the peak with momentum (— 1,0) (degenerate at (|,0)), 
and the peak with momentum (tt, tt) is a shadow of the peak with momentum (0, 0). Haas 
and Dagotto reported the presence of shadow peaks in t-J ladder systems and concluded 
that short-range antiferromagnetic correlations are responsible for the shadow peaks in the 
spin ladder system. 

In order to investigate which of the two parameters U and J is more responsible for 
yielding the shadow peaks, we first set J = and calculate spectral functions for several 
cases of U as shown in Fig. []. We note that the stronger the on-site Coulomb repulsion, the 
broader the bonding orbital becomes, which is more evident for the orbitals, particularly at 
lower values of k x (substantially below the Fermi surface), for example, at k — (0,0) and 
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(f,0) as shown in Fig. |. The overall structure of the bonding and the antibonding bands 
does not change much with increasing U. We find that the shadow peaks can be generated 
even with a small Coulomb repulsion U = 0.1 as shown in Fig. 0(b). While Haas and 
Dagotto's results |TD| are based on the strong coupling limit U ^> t due to the use of the 
t-J Hamiltonian, we discover the shadow peaks in both the weak and the strong coupling 
limits. Second order hopping processes generates an effective interaction of strength ~ t 2 /U 
between the electrons on nearest-neighbor sites. Since the hopping integral along the rung 
(t±) and that along the chain (t\\) are the same, it is hard to determine from their study 
which direction of electron correlations is more important for the shadow peaks. 

In order to thoroughly verify the validity of shadow peaks caused essentially by the 
antiferromagnetic correlation between the two electrons in the rung, namely, by the spin 
singlet state along the rung direction, we computed the spectral functions for the case of 
U = V = 0, Jj_ = J^0 and J\\ = as shown in Fig. f|. We find that the shadow 
peaks appear above the Fermi surface with the Heisenberg interaction along rungs alone as 
displayed in Fig. |](b). Thus the singlet bonding on the rungs is responsible for the presence 
of the shadow peaks. 

Now we reexamine the 5*0(5) symmetric cases of J = 4U. The spectral functions for 
different values of U and J are shown in Figs. |2|(c),(d),(e). As J increases, the energies 
of the bonding and the antibonding orbitals decrease with increasing energy gap. It is to 
be reminded that due to the 5*0(5) constraint to define the relation between J and U, J 
increases more rapidly (4 times faster) than U. Thus it is of great interest to see how the 
spectral positions and energy gaps vary as a result of larger contribution by J. As can be 
seen from the comparison of Figs. |2| and |3], we note that the energy gap is largely contributed 
by J, but not substantially so by U. The dominance of spin singlet states on the rungs for 
the determination of the energy gap is manifest from this study. Indeed this observation is in 
good accordance with the mean field prediction; the energy gap is caused by the Heisenberg 
interaction J as shown in Eq. (|3[). 

The dispersion of the bonding band does not greatly change with increasing values of J, 
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in contrast to that of the antibonding band. As J increases, the antibonding band is lowered 
more rapidly than the bonding band and finally becomes dispersionless (no dependence on 
k x ) as shown in Fig. 0(e). The Heisenberg interaction J, which is taken to be anisotropic 
due to J± 7^ and J\\ = 0, obviously favors the formation of the local rung singlet. The rung 
singlets are clearly local because there exists no phase coherence between the singlets on 
the adjacent rungs, which, in turn, causes the appearance of the dispersionless antibonding 
band. 

We have investigated quasiparticle excitations in the SO (5) symmetric two-leg ladder 
system recently proposed by Scalapino et al. |J The quasiparticle energy dispersion is dis- 
cussed in both the mean field approach and the exact diagonalization study. The dispersion 
consists of two branches, the bonding band with Fermi momentum kp ~ ^ and the anti- 
bonding band with Fermi momentum kp ~ |. The hole spectral functions calculated from 
the Lanczos exact diagonalization method were found to agree well with the mean field 
results. The formation of rung singlet states is responsible for the energy gap in the quasi- 
particle excitations. Finally we find that the presence of shadow peaks above the Fermi 
surface can essentially occur as a result of the singlet bonding on the rungs, indicating 
that the on-site Coulomb repulsion along chains or rungs may not be a prime cause for the 
formation of shadow peaks. 

One of us (S.H.S.S) acknowledges the financial supports of Korean Ministry of Education 
(BSRI-97) and of the Center for Molecular Sciences at KAIST. 
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FIGURE CAPTIONS 



FIG. 1. The quasiparticle energy dispersion obtained by using the mean field approach in the 
weak coupling limit. The parameters are chosen from our Lanczos calculations on a 
2x6 ladder with two doped holes as t\\ = t± = 1, J = 0.4, \ = 0.47, A = -0.15, 
pi = —0.099. B denotes a bonding (k y = 0) band and A, an antibonding (k y = it) 
band. The dashed parts of the dispersions denote shadow bands. 

FIG. 2. Exact diagonalization results of the hole spectral function on a 2 x 6 ladder with two 
doped holes for t\\ — t± — 1, V — 0, and several values of U and J. The dashed 
horizontal line denotes the Fermi energy, and the dotted arrows indicate the shadow 
peaks. The 5 functions have been given a finite width of e = 0.1. 

FIG. 3. Same as in Fig. 2 but for J = 0. 

FIG. 4. Same as in Fig. 2 but for U = 0. 
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FIGURES 
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